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. . Abstract 

"^ ' We study direct energy transfer by multipolar or exchange interactions between diffusive excited donor 

"t^ , and diffusive unexcited acceptors. Extending over the case of long-range transfer of an excitation energy a 
non-perturbative approach by Bray and Blythe [Phys. Rev. Lett. 89, 150601 (2002)], originally developed 

'"O ' for contact diffusion-controlled reactions, we determine exactly long-time asymptotics of the donor decay 

Q ' function in one-dimensional systems. 
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I. INTRODUCTION 

Long-range non-radiative transfer of an excitation energy from excited donor molecules to ac- 
ceptors of the excitation energy is a dominant reaction mechanism in various chemical, physical 
and biological processes [l|, l3) Id) Ij; l5| • To name but a few, we mention fluorescence, luminescence 
or phosphorescence quenching, decay of trapped electrons in glassy media in presence of scav- 
engers, or light harvesting by antennae chlorophyll-b molecules and donation of singlet energy to 
the chlorophyll-a reaction centers in photosynthetic organisms. 

The idea of direct non-radiative transfer has been put forward in the pioneering works of Forster 
and Dexter [2], who determined decay of an immobile excited donor due to dipole-dipole interac- 
tions with immobile, randomly placed acceptors in rigid three-dimensional solutions. Subsequent 
analysis (see Refs.[8| and references therein) extended the consideration of Forster and Dexter 
to arbitrary Euclidean dimensions d and to general forms of donor- acceptor interactions, such as 
isotropic multipolar interactions, for which the rate W{r) of energy transfer is given by 

Wir) = am(^y, (1) 

or interactions mediated by exchange, for which one has 

W{r) = Oe exp (— 7r) , (2) 

where r is the distance separating a given donor-acceptor pair, the constants n, tq and 7 determine 
the interaction type and range (e.g., n = 6 for dipolar, n = 10 for quadrupolar interactions). 

Note that an exponential form in Eq.(l2|) emerges, as well, in another important area - outer 
sphere electron transfer reactions. Kinetics of such electron tunneling processes taking place in 
liquids or glassy media have been also widely studied giving rise to a very beneficial crossfertilization 
of ideas and approaches. 



a 



For the transfer rates in Eqs.([T|) and ^, it was found [8|] that the probability P{t) that the 
donor is still in excited state up to time t obeys, at sufficiently long times, 

P{t) ~ exp [-VdVil - d/n)nArt {cxmtf''] (3) 

for multipolar and 

P{t) ~ exp [-Fd7" V In^ (aet)] (4) 

for exchange-mediated transfer, respectively. In Eqs.(l3|) and ([!]), ua denotes mean density of 
acceptor molecules, Vd = 7r'^/^/r(l -|- d/2) and T{x) is a Gamma-function. The decay forms in 



Eqs.^ 



porous 



and 0) have been also generalized for certain types of restricted geometries - fractals [9 1, 



m 



IC 



n' 



111 ] and various microheterogeneous [l^] media as well as polymer solutions 13l.ll4i|. 
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However, in many situations the donors and acceptors are not immobile. In liquids, both 
donor and acceptor molecules perform diffusive motion. In solids, excitations become delocalized 
because of incoherent hopping between donor sites, which ultimately results in a diffusive transport, 
although the decay kinetics may be still different from that predicted for conventional diffusive 
motion - there always exists a finite probability that an excitation remains on an intitially excited 
donor [l5|, ll6[ . 



It was recognized [13] that random migration of donor and acceptor molecules leads to a 
much more efficient deactivation than the direct transfer between the immobile species; in three- 
dimensions, in particular, one finds that P{t) obeys [171] 



P{t) ~ exp [-47r (Da + Dd) RcsnAt] 



(5) 



where Do and Da are donor and acceptor diffusion coefficients and i?eflf is the effective reaction 
radius. 

Note that the result in Eq.Q has been first obtained for W{p) in Eq.([T]) with n = 6 and 
Da = in Ref . 181] , which analyzed the relaxation of the nuclear magnetization in the presence of 
paramagnetic impurities. Recently, the result in Eq.Q has been generalized to 3D systems with 
donors and acceptors performing anomalous, "fractional" diffusion J19l ]. 

To the best of our knowledge, there do not exist analogous results for low-dimensional, (i.e. ID 
and 2D) systems, although the extension of the approach of Refs.|l7l] for ID and 2D situations 
seems to be straightforward. In doing so, one will find 



Pit) ~ exp [-nA0(t)] 



(6) 



where (/)(t) ~ {{Da + DdW^"^ in ID and (/>(t) ~ {Da + DD)t/ln{{DA + DD)t) in 2D. 

On the other hand, it is clear that neither Eq.([5]) nor Eq.([6]) is an exact solution for the long- 
range energy transfer problem involving diffusive donor and acceptor molecules, but a result of 
a certain assumption. While particles' diffusion coefficients will indeed appear only in the form 
of a sum Da + ^d in the solution of a problem with a single donor and a single acceptor, it is 
not the case in the general situation with a concentration of acceptor molecules. In particular, 
setting Da = 0, one should obtain a crossover to _a_singular behavior characterized by a stretched- 

_^d/{d+2) 



exponential time dependence InP(t) 
in the limit Da -^ 0. 



201 ] , whereas Eq. ([5]) does not show any singularity 



Apart of this, it was recently discovered 2l[ that, remarkably, for diffusion-controlled contact 



C + B ^ B reactions taking place in low dimensional systems (or, generally, in dimension d < df, 



where df is fractal dimension of particles' trajectories in case of subdiffusive motion 23|]), the 
long-time asymptotical form of P{t) is independent of the C particle diffusion coefficient. 



In this paper we extend the non-perturbative approach of Bray and Blythe [2l[ developed 
originally for contact diffusion-controlled reactions to systems with diffusive donor and acceptors 
interacting via distance-dependent isotropic multipolar or exchange transfer rates in Eqs.([T]) and 
([2]). We define, in form of convergent in the limit i — > oo upper and lower bounds, the exact form 
of the excitation survival probability P{t) in one-dimensional systems. More specifically, we show 
that in ID systems, both for multipolar and exchange-mediated transfer P{t) obeys 

fj_\ In P{t) .^ „ / 1 



This exact result proves that, in contrast to predictions based on standard considerations, Eqs.([6]), 
the donor decay function in ID is independent of donor's diffusion coefficient in the asymptotic 
regime. Remarkably, the decay forms appear to be exactly the same as for contact diffusion- 
controlled trapping reactions, (such that reaction takes place upon encounters between particles). 



2l[ despite the fact that in our case reaction proceeds via long-range transfer rates in Eqs.([T]) and 
([2]) and encounters between par ticles do not lead to any particular reaction event. 

In a separate publication [241] , we proceed to show that this is also true for two-dimensional sys- 
tems, while in 3D one may obtain a fluctuation-induced lower bound on the decay function which, 
in some range of parameters, is better (higher) than predictions based on standard Smoluchowski 
approach, Eq.([5]). 

This paper is outlined as follows. In Section II we define the model and introduce basic notations 
for the general d-dimensional case. In Section III we derive a general upper bound on the global 
decay function, while section IV presents the derivation of the lower bound. Next, in Sections 
V and VI, focusing on a ID case, we evaluate the bounds on the global decay functions eplicitly 
for exchange-mediated and multipolar transfer, respectively, and demonstrate that they coincide 
in the asymptotic limit t ^ oo defining in such a way an asymptotically exact result. Finally, in 
Section VII we conclude with a brief recapitulation of our results and an outlook for future work. 



II. MODEL AND BASIC EQUATIONS 

Consider a d-dimensional spherical volume V containing a single excited donor molecule, which 
is initially located at the origin, and K acceptor molecules, placed at random positions. Suppose 
that both donor and acceptors perform conventional diffusive motion with diffusion coefficients Dd 
and Da, respectively. Let the instantaneous positions of the donor and of the acceptors be denoted 
by the (d-dimensional) vectors r(t) and Rfc(t), k = 1,2, . . . ,K. 

We will neglect here the backtransfer to the donor. This neglect is well justified if the donor- 
acceptor energy difference is much larger than ksT, T being temperature and ks the Boltzmann 
constant. We also disregard here donor-specific radiative and radiationless processes. These decay 
channels are independent of the direct energy transfer and thus the overall donor decay function 
factorizes into the product of the donor-specific decay law, exp(— t/r/j), where rpt is the rate of the 
donor-specific decay, times the acceptor determined decay function. Thus we focus here only on 
the non-radiative donor-acceptor transfer. 

One assumes the acceptors to act independently, which means that they contribute multiplica- 
tively to the decay. This assumption is well fulfilled when the density of acceptors is low. Under 
such an assumption, the probability that the donor is still in excited state at time t, for a given 
realization of its trajectory r(t) and given realizations of acceptors' trajectories {Rfc(t)}, is given 
by 

K r ft 



W{pk{t'))dt' 





(8) 



P(r(i),{Rfe(t)}) = J]exp 
fc=i 

where pk denotes the separation distance between the donor and k-th acceptor. 

Experimentally measured property is the global decay function averaged over all possible donor 

and acceptor trajectories 

K 



Pit) = E^ 



n<(o){exp[- r 



W{pk{t'))dt' 



(9) 

R-fe(o) , 



where the symbol E'g {...} denotes averaging with respect to all possible donor's trajectories r(t); 
symbols E^ ,q^ {...} denote averaging with respect to the trajectories of the k-th acceptor, com- 
mencing its motion at position Rfc(O), and finally, the angle brackets stand for the averaging with 
respect to the distribution of the starting positions. Note that presenting P{t) in the form as in 
Eq.([9j), we have already implicitly assumed that all acceptors move independently of each other, 
which is again well-justified for sufficiently low acceptor concentrations. 



After some straightforward calculations, we arrive at the following thermodynaniic-limit expres- 
sion: 



P{t) = E^ {exp [ - nAQ{v{t)- 1)] } 



(10) 



where ua is the mean concentration of acceptor molecules, {ua = K/V when both K,V ^ oo), 
while Q{r(t);t) is the following functional of a given donor trajectory r(t): 

ft 



Q(r(t);t)= / dR(0)ii;^(o) <j 1 - exp 



dt'W{\r{t')-K{t')\) 



(11) 



In the latter equation, R(t) denotes a given trajectory of a single acceptor molecule and E^^q^ {...} 
denotes averaging over all possible trajectories R(t). Note that straightforward averaging in 
Eqs. (|10p and (jlip is a non-tractable mathematical problem since averaging over acceptor tra- 
jectories in Eq. ijlip has to be taken first for a given realization of donor's trajectory and only after 
doing it, one may perform averaging of the exponential in Eq. ijlOp . Consequently, a recourse has 
to be made to approximations. 



III. UPPER BOUND ON THE GLOBAL DECAY FUNCTION: PASCAL PRINCIPLE 



A convenient for our purposes upper bound on the global decay function stems from the so-called 
Pascal principle, which in our terms can be formulated as follows: an excitation on an immobile 
donor molecule survives longer than on a randomly moving one. In other words, P{t) in Eq. ljlOp 
is bounded by 



Pit) < Pu{t), 



(12) 



where Pu{t) describes the decay of an immobile donor, fixed at the origin, due to a concentration 
UA of diffusive acceptor molecules, 

rt 



Pu{t) = exp 



-HA I dK{0)E^^,^ 



1 — exp 



dt'W{\K{t')\) 



(13) 



The inequality^ 
proven in Ref. 



n 



in Eq. (jl2p has been first conjectured in Ref.[21l| for contact trapping reactions and 
2j] for one-dimensional systems. In Ref. 25|], Eq. (jl2p has been proven for a rather 



general class of random walks on d-dimensional lattices. We also remark that a similar statement 
has been proven earlier in Ref. [l5|] for the process of an excitation energy migration via distance- 
dependent transfer rates on a disordered array of immobile donor molecules and quenched by 
randomly placed immobile acceptors. It was shown that the survival probability of an excitation 



can be only decreased because of random motion not correlated with acceptors' spatial distribution. 
However, no rigorous proof of such a statement exists at present for diffusion-controlled long-range 
reactions although it is intuitively clear that the inequality in Eq. (jl2p should hold in this case too. 
We thus assume, without proof, that the inequality in Eq.([T2]) is also valid for the model under 
study. 



Next, applying Feynmann-Kac theorem 27|, 



28t | one may show that 



^R(O) \ exp 



dt'Wi\Kit')\) 



dRGi(R|R(0)), 



(14) 



Gt(R|R(0)) being the Green's function solution of the following Schrodinger equation: 
d 



dt 



Gi(R|R(0)) = Z)AARGi(R|R(0)) - VF(|R|)Gi(R|R(0)), 



Gi=o(R|R(0)) = <5(R-R(0)), 



(15) 



where Ar, is a d-dimensional Laplace operator. 

Note that Eqs. ()15p presumes that donor and acceptors are point-like, non- interacting particles. 
In reality, they possess hard-cores and can not approach each other at distance less than a, equal 
to the sum of donor and acceptor radii. This means that Eqs. pSj) are to be complemented by a 
reflective boundary condition at |R| = a [4]. 

Taking advantage of Eqs. (|14p and (jlSp . we can formally rewrite Eq. (jl3p as 



Pu{t) = exp 



UA I dt' I (IK * ^ ^ 



Gt(R) 



dR(0)Gt(R|R(0)). 



(16) 



/o J dt' 

Assuming next that Gt(R) is independent of angular variables such that G'j(R) = Gt{r), where 
r = |R|, we get the following compact expression: 



Puit) = exp -UA / dt'kuit') 
Jo 

in which equation ku{t) is determined by 

fOO 

ku{t) = dVd r^-^W{r)Gtir), 

J a 

and Gt{r) obeys 



(17) 



(18) 



dGtjr) 
dt 

Gt=o 



j^^. d^Gtjr) ^ d-ldGt{r) 
"^ "^ r dr 

dGtir] 



Qj.2 

1; Gi(r^oo) = l, 



dr 



W{r)Gt{r) 
= 0. 



(19) 



Equations (fTTj) . (fTHj) and (fT9]) thus define the upper bound on the global decay function Pit) in 
systems with diffusive donor and acceptors. 



IV. LOWER BOUND ON THE GLOBAL DECAY FUNCTION. 



We turn now to the derivation of a lower bound on P{t) in Eq. (|lUp . Following Ref. 21[ (see also 



Ref. 



26t|). we make the following steps: 



(i) suppose that for a given initial placement of acceptors, a closest to the origin acceptor appears 
at distance I. Thus, a notional spherical volume Vi of radius /, centered at the origin, is initially 
completely devoid of acceptors. 

(ii) performing averaging over donor's trajectories {r(t)}, we consider only such trajectories which 
never leave Vi up to time moment t. Since Q{r(t);t) in Eq. (jlip is always positive definite for any 
particular realization r(t), such a constraint naturally leads to a lower bound on P{t), i.e. 

^^ {exp [ - nAQ{r{t); t)] } > ^o^r(i)eV, {exp [ - nAQ(r(t); t)] } , (20) 

where E^^^^^^y {...} denotes averaging over a subset of all possible donor's trajectories such that 
they do not leave Vi during time t. 

(iii) considering the term responsible for long-range transfer, Q{r{t);t), we suppose that the donor 
is always located on the surface of Vi at position closest to the instantaneous position of the 
acceptor. Since W{p) is a strictly decreasing function of p, for any r{t) E V/, one has H^(|R(t)| — /) > 
VF(|r(t) — R(t)l) and hence, Q{r{t);t) can be majorized by 



Q(r(t);t)<Q(/;t) 



dK{0)E^, 



(0) 



1 — exp 



dt'W{\K{t')\ -I) 



(21) 



Note now that the right-hand-side of the inequality in Eq. (|2ip is independent of the donor's 
trajectories. 

Consequently, collecting (i) to (iii), we arrive at the following lower bound on the global decay 
function 



P{t) > Pvoid(0 X ^o'^rWev, {1} X exp [-nAii(/;t)] 



(22) 



In this equation Pymdi^) is the probability of having an acceptor-free spherical void of radius I. For 
random intial placement of acceptors, one has 

id 



Pvoid{l) ~ exp 



UAVd r 



(23) 



Further on, in Eq. (j22p the symbol -E^r (,\^y {1} denotes the measure of such donor's trajecto- 
ries, which commence at the origin and never leave V/ during time t; at sufficiently large times, 
{1} is given by 






^o,r(i)ev, {l}~exp 



-Zd- 



,DDt 



p 



(24) 



Zd being the first zero of the Bessel function J[d-2)/2{^)- 

Combining the expressions in Eqs. (|23p and (|24p . and assuming spherical symmetry, we finahy 
obtain 



In the latter equation, 



-UAVdl" 



z^^ UA / dt'ki{t') 



P 



h{t)=dVd r^-'W{r-l)Gt{r)dr, 

Jl+a 



(25) 



(26) 



while Gt{r) is the solution of 

dGtjr) ^ 
dt 



Da 



d^Gtir) d-ldGt{r) 



Qy.2 



+ 



dr 



W{r-l)Gt{r), 



(27) 



Gt=o{r) = 1; Gt{r ^ oo) = 1, 
subject, in virtue of condition (iii), to a reflection boundary condition imposed at r = / + o: 

= 0. (28) 



dGAr) 



dr 



r=l+a 



Equations (|25p .(j26 p . (|27p and (|28p define a family of lower bounds on the global decay function in 
systems with diffusive donor and acceptors, dependent on the radius / of the notional volume Vi 
encircling the donor and devoid of acceptors. 

To get the optimal lower bound, we will have, in the usual fashion, to maximize the result 
with respect to /. Below we consider lower and upper bounds on the global decay function in 
one-dimensional systems with long-range transfer (Eqs.([T]) and ([2])) between diffusive donor and 
diffusive acceptors. Corresponding results for two- and three-dimensional systems will be presented 
elsewhere J24|. 



V. ONE-DIMENSIONAL SYSTEMS: EXCHANGE-MEDIATED TRANSFER. 



A. Upper bound. 



Consider first the derivation of an upper hound in one-dimensional systems with a transfer 
mediated by exchange. Here, Laplace-transformed with respect to time variable t, solution of 
Eqs.dSZ]) and ^ reads: 



Gx{r) 



r-oo 

/ dtexp[-At]Gt(r) = Cil^{x) + C2K^{x) + 
Jo 



+ 



xr(i-f) /x^f '■' 



2 



(xOe 



1-e 



2^f-l 



d^, 



(29) 



where Ky.(x) and Ix{x) are modified Bessel functions, the integral term in the second-hne is a 
particular solution (Lommel function) and 



r, , a, 2 / Og , 2 / A 



x = a;exp[-7-], xq = u;exp[--f-], uj = -J — , and x=-yjy- 



(30) 



Now, note that as r ^ oo, x ^ 0, I^i^) -^ 0, the last term on the rhs of Eq. (129p tends to 1/A, while 
Ky.{x) diverges. Hence, we set C2 = 0. Further, we get that the reflective boundary condition at 
the closest approach distance is fulfilled when 

2xr(i - f ) 



Ci 



^Y' l"h^,ixoc)e-'Hi-e)"-'dc (31) 



A(/x-i(^o) + ^x+i(^o)) V2y Jo '2 
Plugging Eqs. (i29]) and (i3T]) into Eq. ifTSJ) and performing intergation, we find that the Laplace- 
transformed reaction constant ku{X) is given by 



kuiX) 



DaIxI 



A 



2F3I 1,1:2,1 



X ^^X xl 



2' 2' 47 (i + |)(i_2d)r(l + ;^) 



XqY+X 



1+. 1F2 2;2-f,2 + 



X. £0 

2' 4 



1F2 l + f;2 + f,l + x;^ 



^x-i (2^0) + -^x+i (2^0) 
where pFg denote generalized hypergeometric functions. 

Leading small-A (large-t) asymptotic behavior of fcu(A) in Eq. (f32]) follows 

1 



(32) 



ku{X) ~ 2 



Da 

A I + VtTA' 



(33) 



where 



Ki (xo) + (1/2 - g + In (2/xo)) Ji (xp) 
^1 (2;o) 

C ~ 0.577 being the Euler constant. 

This yields, in t-domain, the following asymptotical behavior 



I^Da' 



(34) 



dt'ku{t') 



4,/^ I 1 



vr 



4t V t 



(35) 



Consequently, in Id systems with transfer mediated by exchange we have the following upper bound 
on the global decay function: 



P{t) < exp 



.^nA^ + 2nA^/D^e + o(^-^^ 



(36) 



Before we proceed to the derivation of the lower bound, a few comments are in order: 

(a) first of all, we notice that the right-hand-side of Eq. ()36p coincides with the solution of the 



10 



so-called target problem - probability that an immobile target survives, in one dimension, up to 
time t in presence of diffusive scavengers which may "destroy" the target upon the first encounter 



with it 29|. Therefore, in one dimension, at sufficiently long times the kinetic behavior of long- 
range transfer proceeds exactly in the same way as for contact diffusion-limited target annihilation 
reaction, despite the fact that here the boundary condition imposed on the donor's surface is 
reflective and the deactivation of the donor happens, at rate aeexp(— 7r), at any donor-acceptor 
distance r. 

(b) parameter Tg in Eq. (j34p is the crossover time to the asymptotic stage In Pu{t) ~ — t^'^ for 
exchange-mediated transfer in one-dimensional systems with immobile donor and mobile acceptors. 
Note that ^'^DaTs is a non-monotonic function of xq. It is large ~ l/a^o when xq <^ 1 (i.e., when 
Da is large), such that Tg ~ ^'^DA/a'^- In this case, one would first observe, for < t < Tg, 
an intermediate asymptotical behavior In P(t) ~ —jat, which will then cross to the asymptotical 
behavior in Eq. (j36p . Next, note that ^"^DaTq is also large when xq ^ 1, which happens when 
Da is small. Here, Tg ~ In {xq/2)/j'^Da, i.e. Tg is proportional to the first inverse power of Da 
(with logarithmic corrections). In this case, the asymptotic decay in Eq. ()36p succeeds the static 
quenching decay in Eq.Q, which is valid in progressively larger time domain the closer Da is to 
zero. 

(c) finally, we remark that despite the fact that the result in Eq. (|35p is independent of both ag and 
7, which are the only parameters characterizing the transfer rate and thus "represent" reaction, it 
does not mean that it can be simply obtained by expanding G{r) = Y2'n=o '^"^'^(r) and considering 
the zeroth term only. In general, Eq. (l35p is essentially a non-perturbative result and can not be 
obtained using a perturbative expansion of Gt{r) in powers of ag, unless, of course, one manages to 
sum the whole series. On the other hand, Eq. (j35p can be straightforwardly derived approximating 
the transfer rate by a step- function ("square well" approximation). 

Indeed, suppose that 7a < 1 and consider separately solution of Eqs. ([27|) and (p8]) for a < 
T ^ 1/7 and r > I/7. In the first interval we approximate exp[— 7r] by exp[— 7a], and find that 
the Laplace-transformed solution of the Schrodinger equation which obeys the reflecting boundary 
condition reads 



.(1) ^ L 

A -|- agexp[— 7a] 



G\^ = ^ , ^^ ^ + Ci cosh 



A -|- Ogexpi— 70] 

— [r — a) 



(37) 



Da 

On the other hand, in the domain r > I/7, the transfer term can be neglected, and we have 

(38) 



A 



GrW = T + ^2exp 



11 



Da\ 1 



Since Gx{r) and its first derivative have to be continuous functions at r = I/7, we have two 
complementary equations which define the coefficients Ci and C2. Determining these coefficients, 
we find that the leading small- A behavior of ku{X) follows 



/•1/7 
A;„(A) = 2aeexp[-7a] / G 

J a 



!"■ 



r)dr 



Da 



1 



A i + yip 



where 



ri = ^^^P[^coth2(f 



Op 



(39) 



(40) 



Note that A;„(A) in Eq. (j39p has exactly the same form as ku{X) in Eq. (j33p . which means that the 
"square well" approximation captures well the leading behavior of the effective reaction rate. The 
crossover time Tg has a different form compared to the exact one, Eq. ()34p : it exhibits, however, 
quite a "correct" behavior in the case xq <C 1 (fast diffusion) when Tg ~ 7^D^/ag. 



B. Lower bound. 



Consider now a lower bound on -P(i) for one-dimensional systems with transfer mediated by 
exchange interactions. Laplace-transformed solution of Eqs. (|27p and (j28p reads 



<;.M = c,/,M) + ^^I(l^(=i)'^^ 



2 



{x0)e' 



Hl-e)^''d^, (41) 



where Ci is given by Eq. (l3ip . Plugging the expression in Eq. ()4ip into Eq. ()26p and performing 
integration, we find that ki{X) obeys 



k{X) = A:„(A), 



(42) 



where A;„(A) is determined by Eq. (j32p . Consequently, the lower bound on P{t), Eq. (j25p . at suffi- 
ciently long times attains the following form: 



Pi{t) ~exp 



-2nAl-.^£^-,nJ^-' 



P 



IT 



(43) 



As we have already mentioned, the result in Eq. ()43p represent rather a family of lower bounds 
dependent on parameter / - radius of a notional volume initially devoid of acceptors. The "best" 
lower bound thus would be the highest one. Optimizing Eq. ()43p with respect to /, we find that the 
highest lower bound is achieved when / = (Tr^D^it/n/i)^''^, and is given by 



Pi,max{t) ~ exp 



-AnAj^-3nf{.'Dntf' 



TT 



(44) 
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On comparing the asymptotic behavior predicted by the maximal lower bound in Eq. (j44p against 
the upper bound in Eg . ([55]) we notice that both bounds converge asymptotically to give an exact 
result in Eq.©. 



VI. ONE-DIMENSIONAL SYSTEMS: MULTIPOLAR TRANSFER. 



A. Upper bound. 



Consider now, within the "square well" approximation, an upper bound in case of multipolar 
transfer in Eq. ([!]). Approximating the actual transfer rate W{r) in Eq.([T]) by a step-function 



W{r) 



am{ro/a)'^,a < r < ro, 
0, r > ro 



we find that in the interval a < r < tq Laplace-transformed solution of Eq. (|19p obeys 



G 



(1) 



1 



while in the domain r > ro it follows 



+ Ci cosh 



\ + am{ro/a) 



Da 



-(r-a) 



(45) 



Gf(r) = l + C2exp 






(46) 



Constants Ci and C2 are to be chosen in such a way that both Gx{r) and its first derivative are 
continuous functions at r = vq. 

Determining these constants, plugging Eq. (j45p into Eq. (jl8p and performing integration, we find 
that the Laplace-transformed A;u(t) is given by 

-1 



'"W = ^>/^l'^,„.... 



Aa» ,oth(./ °"''y°'" (ro-a) 



m' 



Da 



2a^(ro/a)"(ro - a) 

A + a™(ro/a)" ^ ' 



This yields, in t-domain 

/ K{t')dt' 
Jo 



4./^ 11 



vra" 



vr 



AamTQ t 



■ coth 



'^^,.-,).o(l 



(48) 



Consequently, the global decay function P{t) in one-dimensional systems with diffusive donor and 
acceptors interacting via multipolar transfer rate in Eq.([T|) is bounded from above by 



P{t) < exp 



-AuAxI 1" 2n^4 / coth 



'am(ro/a)" 



vr 



O^m^Q 



Da 



{ro - a) 



(49) 



13 



B. Lower bound. 

Turning next to evaluation of the lower bound on P{t) we introduce parameter 5 > a, and 
approximate the actual transfer rate by a step-function of the form 



W{r) 



am(?"o/a)", l + a<r <l + 5, 
0, r>6 



Approximate solution of Eqs. (l27p in the interval l + a<r<l + 6 has the form 



G 



(1) 



1 



X + am{ro/a) 



+ Ci cosh 



'\ + am{ro/ay 



Da 



{r — I — a) 



(50) 



while in the domain r > / + (5 it is given by 

Gf{r) = \ + C2e^V 



A 



DiC--'-*) 



(51) 



Again, requiring continuity of G\{r) and of its first derivative at r = / + 5, we determine Ci and 
C2, which yields, after straightforward calculations, the following expression 

-1 



^■'<^'-^/^(-\/^ 



coth 



1^^^'-'^))) + A + a.(Wa)" ' ^''^ 



We find then that in the t-domain, the leading behavior of L ki{t')dt' is given by 



„M.M*-.v^-VS:^"M\/'^<^-"'i^''^^ 



(53) 



Consequently, an optimized lower bound on P{t) reads 



P{t) > exp 



4nA\l^-3nf{7:'Dntf' + 2nA. 



fDAa"^ 



vr 



n 
ml 



coth 



'am(?^o/a)" 



CXmTi 



Da 



{S-a) 



(54) 



On comparing the results in Eqs. (|49p and (j54p . we notice that again both bounds converge as 
t — > 00 determining exact asymptotic decay of the excited donor, Eq.©. 

VII. CONCLUSIONS 

To conclude, we have studied analytically direct energy transfer between diffusive excited donor 
and diffusive unexcited acceptors mediated by multipolar or exchange interactions. Extending a 
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non-perturbative approach by Bray and Blythe [21] (originally developed for contact diffusion- 
controlled reactions) over the case of long-range transfer, we have determined exactly long-time 
asymptotics of the donor decay function in one-dimensional systems. We have shown that the 
leading long-time behavior is independent of the diffusion constant Do of the donor molecule, and 
has exactly the same form as that descibing contact process. This finding is in aparent contradiction 
with the results in Eqs.([5]) and ([( 



We proceed to show elsewhere |24| ] that also in two-dimensional systems the leading long-time 
behavior will be independent of Do, while in 3D a similar approach will give rise to a fluctuation- 
induced lower bound on the decay function which, in some range of parameters, is better (higher) 
than predictions based on standard Smoluchowski approach. 
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